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Nonleptonic decays of charmed mesons into two pseudoscalar mesons or one pseudoscalar meson 
and one vector meson are studied on the basis of a generalized factorization method considering the 
resonance effects in the pole model for the annihilation contributions. Large strong phases between 
different topological diagrams are considered in this work, simply taking the phase in the coefficients 
Oi. We find that the annihilation- type contributions calculated in the pole model are large in both of 
the PP and PV modes, which make our numerical results agree with the experimental data better 
than those previous calculations. 
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I. INTRODUCTION 



Nonleptonic decays of charmed mesons are very interesting as they can provide useful information on flavor mixing, 
CP violation, and strong interactions [l|. They may also shed light on any new physics signal through D° — D° mixing 
and rare decays [^-Q]. The CLEO-c and two B factories experiments have given many new results on this subject and 
more are expected soon from the BES-III experiment. Besides, theoretical studies have been in progress for decades. 

With the heavy quark effective theory, many QCD-inspired approaches, such as the QCD factorization approach 
Q , the perturbative QCD approach Q , and the soft-collinear effective theory [l(j , successfully describe the hadronic 
B decays. However, this is not the case for the D meson decays. These approaches do not work well here, due to the 
mass of charm quark, of order 1.5GeV, which is not heavy enough for a sensible heavy quark expansion, neither light 
enough to apply the chiral perturbation theory 

After decades of studies, the factorization approach is still one of the effective ways to deal with the two-body 
charmed meson decays[llj. However, it is well known that some difficulties exist in the naive factorization approach: 
the Wilson coefficients a±(fj,) and 02 (/i) of effective operators are renormalization scale and 75-scheme dependent; 
and the color-suppressed processes are not calculated well due to the smallness of 02, etc. In order to solve these 
problems, the so-called generalized factorization approach was proposed [l2|. The Wilson coefficients ai(fi) and a<z (/z) 
are not from direct calculation any more but are effective coefficients to accommodate the important nonfactorizable 
corrections fl3l ]. In these naive or generalized factorization approaches, there is almost no strong phase. But large 
strong phases have been found from D decay experiments. A corresponding large relative strong phase between the 
factorization coefficients a\ and 02 has been discussed in flil Il5|. 

Unsatisfied by the short comings of the factorization approach, the model- independent diagrammatic approach, with 
various topological amplitudes extracted from the data, is recently applied to two-body hadronic D decays |16l - [l9| . 
They use SU(3) symmetry in their analysis to avoid model calculations. All the parameters are fitted from experiments 
to give a better agreement with the experimental data but with less predictive power. More precise predictions are 
limited in this approach due to the uncontrolled SU(3) breaking effect These analyses also show that large 
annihilation type contributions are needed to explain the data, which can not be calculated in the naive or generalized 
factorization approaches. The kind of pure annihilation type D meson decays also needs to be systematically analyzed 

M- 

In another aspect, the hadronic picture description of nonleptonic weak decays has a longer history, because of their 
nonperturbative feature. Based on the idea of vector dominance, which is discussed on strange particle decays 2l|, the 
pole-dominance model of two-body nonleptonic decays is proposed [22[ ■ Beyond the vector dominance pole, this model 
also involves scalar, pseudoscalar, and axial-vector poles. For simplicity, only the lowest-lying poles are considered. 
This model has already been applied to charmed meson and bottom meson decays [22|, [23| , where it is approved that 
this model is more or less equivalent to the factorization approach in the first order approximation. 

In this work, we will use the generalized factorization approach but with a relative strong phase between the Wilson 
coefficients ai and 02 to accommodate the nonfactorizable contributions. For the uncalculablc annihilation type 
contributions, we use the pole model. In this case, we can really calculate most of the important contributions in the 
hadronic D decays, which are demonstrated in the model-independent diagrammatic analysis [Tg| . 

The outline of this paper is as follows. In Sec. II, we give the formulas of this work, showing the generalized factor- 
ization approach and the pole-dominance model. In Sec. Ill, our results are given and compared to the experimental 
data and those of the diagrammatic approach and the calculations considering the final-state interaction of nearby 
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FIG. 1: Emission-type diagrams in the factorization approach. 



resonances effects. Summary and conclusions are followed. 

II. FORMALISM 
A. The factorization approach 

First we begin with the weak effective Hamiltonian H e f f for the AC = 1 transition [24[ : 

n eff = ^VcKM(c 1 o 1 + c 2 o 2 ) + h.c., (i) 

where Vckm is the corresponding Cabibbo-Kobayashi- Maskawa (CKM) matrix elements, Ci,2 are the Wilson coef- 
ficients. The current-current operators are 

01 = u Q 7 M (l - 75)92/} • g3/37 M (! - 75)c Q , 

02 = S Q 7 M (1 - 7 5 )<72 Q • g3/37 M (! ~ 75)c/3, (2) 

where a, /3 are color indices, and 92,3 are d or s quarks. 

The color-favored emission diagram corresponding to D — > PP decays, with P representing a pseudoscalar meson, is 
shown in Fig. [T](a). Under the factorization hypothesis, the transition matrix element of hadronic two-body charmed 
meson decays is factorized into two parts [ll| : 

{P l P 2 \U e{f \D) = ^|Uc/o/ai(P 2 |S7 M (l-75)<72|0)(Pi|g37' i (l-75)c| J D). (3) 

Similarly, the contribution of the color-suppressed diagram shown in Fig. [1] (b) is given as 

(P 1 P 2 \H 6 f f \D) = ^VcKMa2(Pi\q 3 ^(l-J 5 )q2\0)(P2\ur^~^)c\D), (4) 

where 

ai (0) = C a M + ^, a2 (/i) = 1 (/i) + ^, (5) 

correspond to the color-favored tree diagram (T) and the color-suppressed diagram (C) respectively in the naive 
factorization, with the number of colors N c = 3. They are assumed to be universal and process- independent in the 
native factorization approach. The current matrix elements in Eqs. (|3l4[ ) arc evaluated in terms of transition form 
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factors and decay constants. For D — > PP decays, the form factor is denned as follows: 



(P(k)\q 3lfl (l - l5 )c\D(p)) 



(P + *)p 2 % 



q 2 



F^ p (q 2 ) 

e^o^V), (6) 



where q = p — k, and F^ are the corresponding transition form factors. The decay constants fp of pscudoscalar mesons 
are defined as 

(P(?)|? 1 7m(1-7S)«i|0> = ifpq„. (7) 
In terms of decay constant and transition form factors, the decay amplitude of Figs|T] (a) and (b) are then 

{P l P 2 \U eff \D) T = i^VcKMaifPziml-mpJF^impJ, (8) 

(P x P 2 \H ef ;\D) c = i^VcKMa2fpMl-™p 2 )F^ P2 (mpJ. (9) 

The diagrams for D —> PV decays, with V denoting a vector meson, are shown in Figs. [T] (c) and (d). If the 
out-emitted particle is a pscudoscalar meson, the matrix element of Fig. [T] (c) is 

(PV\H eff \D) = ^VcKMa 1 (P\u^O--j 5 )q 2 \0}(V\q 3 ^(l--f 5 )c\D}. (10) 

The D — > V transition form factors are usually defined as 

(V(k)\q 37fl (l- l5 )c\D(p)) = \ e^ pr7 s*»pOk°V(q 2 ) 



q 



+i ( (P + *)„ - ) -^-^A^(q 2 ) 

q A J mpt + my 



where e* is the polarization vector of the vector meson, and A, and V are corresponding transition form factors. 
Utilizing the form factor definitions we get the result for Eq.(fTT 



(PV\H eff \D) = ^V C KMaifpmvA^ v (m 2 p)2(e*-p D ). (12) 

If a vector meson is out emitted, the matrix element of Fig. QJc) is 

(PV\H ef f\D) = ^VcKMai(V\u^(l- l5 )q 2 \0)(P\q 3 ^(l-7 5 )c\D) 

= -^VcKMaifvm v F^ p {m 2 v )2(e*-p D ), (13) 

where the decay constants fy of vector mesons are defined as 

(ng)l?i7M(l-75)g 2 |0) = fymye;(q). (14) 

For the color-suppressed diagram, FigfjTd), we have similar formulas as Eqs. (|12ll3[) . but with the Wilson coefficient 
changed from a\ to a 2 . 

As mentioned in the Introduction, the Wilson coefficients a\ and a 2 are renormalization-scalc dependent in the 
naive factorization approach and it fails to describe the color-suppressed processes with too small a 2 « —0.1. So the 
generalized factorization method is proposed to include the nonfactorizable contributions [l3| , 

,eff 



af = G a 0i) + (Jj- + , of = + C 2 (/x) (1- + X2 (m)) , 
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where the terms Xi characterize the nonfactorizablc corrections involving vertex corrections, hard spectator interac- 
tions, final-state interactions, resonance effects, etc. These Xiil 1 ) wm compensate the scale- and scheme-dependence 
of the Wilson coefficients, so that a,'s are physical now. Without confusion, we will drop the superscript "eff" in the 
effective Wilson coefficients for convenience in the following discussions. In the large- N c approach, the 1/N C terms 
are discarded [25j. equally with a universal nonfactorizable term Xi = X2 = —1/N C , hence, 

01 C 2 (m c ) = 1.274, a 2 w Ci(m c ) = -0.529. (16) 

This implies a null relative strong phase between the two kinds of contributions. However, the experimental data tell 
us that there should be a large strong phase between a\ and a 2 . On the other hand, the existence of relative phases 
is reasonable for the importance of inelastic final-state interactions of the D meson decays, in which the on-shell 
intermediate states contribute imaginary parts. Therefore, we consider a relative phase between the coefficients a\ 
and a 2 in this work, so that 

ai = \ai\, a 2 = \a 2 \e iS , (17) 

where we set a\ real for convenience. 



B. Pole-dominance Model 



The annihilation type diagrams are neglected as an approximation in the factorization model. However, considerable 
contributions come from the weak annihilation diagrams in the D decays, which can be demonstrated by the difference 
of lifetime between D° and D + . Hence, we will calculate them in a single pole-dominance model. For simplicity, only 
the lowest-lying poles arc considered in the single-pole model. Taking D° — > ir + K*~ as an example, the annihilation- 
type diagram in the pole model is shown in Fig. [2ja). D° goes into K° via the weak interaction in Eq. ([I} shown in 
terms of quark lines in Fig. EJb), and then decays into ir + K*~ through the strong interaction. Angular momentum 
should be conserved at the weak vertex and all conservation laws be preserved at the strong vertex. So it is a 
pscudoscalar meson as a resonant state for D — > PV decays. The weak matrix element is evaluated in the vacuum 
insertion approximation [23j . 

(K°\U eff \D°) = ^y; s K lrf ar(K | S - 7M (l- 75 )d|0)(0|^(l-75) C |^ ) 

= ^V; s V ud a^ v f K f D ml. (18) 

where the subscript E of the Wilson coefficient a^ v denotes a W-exchange diagram for D — > PV, otherwise a F A v 
corresponds to W-annihilation contributions. In fact, the effective Wilson coefficients of the W-annihilation diagrams 
and W-cxchange diagrams have the same form as a\ and a 2 in Eq. (|15[) . that is 

a E =Ci(//) + C 2 (m) (Jj- + Xe(v)J , aA = C 2 ((J,) + C 1 ( f i)(J^+XA(^ , (19) 

where Xa(e) represents the nonfactorizable contributions in the annihilation (exchange) process. Since the nonfactor- 
izable contributions in these kinds of diagrams are large and with relatively different strong phases, we use different 
symbols to avoid confusing in our approach for these collective effective Wilson coefficients as cle and a a , respectively. 
Strong phases relative to the emission diagrams are considered in the Wilson coefficients. The effective strong coupling 
constant of K° to ir + K*~ is defined through the Lagrangian 

Cvpp = igvppV^P^^P), (20) 
where gvpp is dimcnsionless. Inserting the propagator of the intermediate K° meson, the decay amplitude is 

(n+K*-\H eff \D°) = ^V*V ud a E v f K f D g K , K « ™* D 2 2(s* -p D ). (21) 

Similarly, in D —> PP decays, it is a scalar meson as a resonant state. The effective strong coupling constant is 
described by 



£spp = -gsppmsSPP, 



(22) 
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FIG. 2: Annihilation diagrams in the pole model. 



where ms is the mass of the scalar meson. Besides, the scalar meson decay constant of the vector current is defined 
as 

(S{p)\q 2 ^qi\0) =fsP M . (23) 
Therefore, the corresponding matrix element is 

(PP\He ff \D) = -i^VcKMa A (a E )f s f D g SP p . (24) 

V2 m z D - nig 

As a convenience of reference, the various decay formulas of individual decay modes are collected in the appendix. 
Note that some of the intermediate resonances are unstable particles which have large width and therefore contribute 
large relative phases. These phases are absorbed in the effective Wilson coefficients a a and oe for convenience. 



III. NUMERICAL ANALYSIS 



A. Input Parameters 



In order to calculate the emission type diagrams in the factorization approach, we need to know the transition form 
factors and meson decay constants. The decay constants of 7r, K, D and D s are taken from the particle data group 
(PDG) (26[, others are from (2?|], all of which are summarized in Table Q] There exist many models to parametrize 
the transition form factors and their q 2 dependence [28l - [46j . In this work we shall use the dipole model [2|| : 



where m po i e is the mass of the pole. The corresponding poles are D* for F^' Dv ,D,K , D* for F^^' " 11 ^ jj ^ 
A Q p,Du '> Ds } and D s for A^ K ' Ds . The transition form factors and parameters of D to it and K arc taken from 
the recent CLEO-c measurement [29[ , D — > r\ q are from [3(j , and others from [28| , all of which are shown in Table HU 
For the final states involving r\ or rf , it is convenient to consider the flavor mixing of r/ q and rj s with a mixing angle 



V \ _ ( cos V ~ sm <P \ I Vq 



rf j \ sm <p cos <p I \ rj s 



(26) 



where r\ q and rj s are defined by 



rj q = -j={uu + dd), rj s = ss. (27) 
v2 



A recent experimental measurement from the KLOE collaboration gives the mixing angle (f> = (40.4 ± 0.6)° j47j|. 
The decay constants of rj or rf arc defined by 

(0|u7 M 75u|7?(p)) = if^Pn, (0\djp,j5d\r){p)) = iffon, (0|s7/*7ss|j7(p)) = */^V> 
(0\u W u\ V '{p)) = if;, P „ (0\d W d\r)'(p)) = if* Pll , (0\s W s\ V '{p)) = if°,p„ (28) 
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TABLE I: Meson decay constants (MeV). Those of ir, K, D, and D s are from PDG[H, others are from [13|. 



fir fK 




/a- 


/« U fD 


fD, 


130 156 


216 


220 


187 215 207 


258 



TABLE II: The D meson transition form factors and dipole model parameters Qi,2- The parameter ot\ = «2 + 1 if only q?2 is 
shown in the table. The form factors and cti parameters of D to it and K are from [2!|, D — >• 77, from [33], and others from 





Po 




P DK 


pDK 


f d s k 


p D sVs 


F(0) 


0.67 


0.67 


0.74 


0.74 


0.72 


0.78 


(22 


0.21 


0.24 


0.30 


0.33 


0.20 


0.23 




A Dp 


aDu> 


A DK* 


4 DsA " 


aDs^ 
A o 




F(0) 


0.66 


0.66 


0.76 


0.67 


0.73 




Q?2 


0.36 


0.36 


0.17 


0.20 


0.10 






^0 


^0 


D s ^ Va( M v , 
-M) 


1 pD^lq 


pD^Vq 




F(0) 


0.72 


0.78 


0.78 


0.69 


0.69 






0.41 


0.33 


0.21 


1.03 


0.39 




OL2 


0.70 


0.38 


0.76 


0.29 


0.01 





where 

1 1 

fU_ fd _ rq ru fd rq /r,Q\ 

J v J n -s/2 ' '' '' \/2 



With the ansatz in 1481 , we have 



/« = / g cos </>, /* = -f a sin 0, 

/«, = /, sin 0, f°,=f s cos<t>. (30) 

It is assumed that / 9iS obtained from the ?y 9iS components of the wave functions are independent of the meson involved. 
We use that f q = (1.07 ± 0.02)/* and f s = (1.34 ± 0.06)/* from [H|. 

The form factors of Z? — > ry g in Table [TT1 denote that of D — >• 7/, ufi d g, not Z? — ^ (uu + dd), hence, 

F°- ,ri = ^F D ^ Vq cos (f> : F Ds ^ r > = -F Ds ^'i° sin0, (31) 
F D ^' = ■1=F D ^« sin0, F Ds ^' = F Ds ^ cos<f>. (32) 

In order to calculate the annihilation-type diagrams in the pole model, we have to know the effective strong coupling 
constants between the intermediate state and two final states. Some of them are obtained directly from experiments. 
Some others are related to the known ones using SU(3) symmetry. Although the intermediate states are a little off 
shell, in the pole model they are used as on-shell resonant states. So these on-shell strong couplings are used to 
calculate the annihilation diagrams in an approximation. 

There are many scalar mesons discovered by the experiments. The existence of the lightest scalar nonet with the 
mass smaller than or close to 1 GeV has been a problem for many years [2^. It is still controversial that they arc 
primarily the four-quark bound states or two-quark scalar states. In this work, we use ^(1430), ao(1450), /o(1370), 
and /o(1500) as intermediate mesons in the pole model for D — > PP decays. The decay constant of if „ is calculated 
in several methods, such as the finite-energy sum rule[4!|, the generalized Nambu-Jona-Lasinio model [5(|, and so on. 
We shall use the results in (49j . 

f K . = (42 ± 2)McV. (33) 

For all other scalar mesons, we take the same value of decay constants as Kq in the flavor SU(3) limit for simplification. 

The corresponding effective strong coupling constant between Kq and the final states of Kir is evaluated by Qk* } k tt = 
2.7 from the decay of A"q(1430)° — ^ ir~K + . Other couplings of gspp are of the same value in the SU(3) limit. For 



7 



TABLE III: Branching ratios for Cabibbo-favored decays of D — > PP(%). The predicted branching ratios with both 
annihilation- and emission-type contributions (all) and with only emission-type contributions (emission) are given together 
with the experimental dataJH^j, the recent results from the diagrammatic approach and the calculations considering the 
final-state interaction (FSI) effects of nearby resonances [5lT[ as comparison. 



Modes 


Br(FSI) Br(diagrammatic) Br(emission) 


dl i an i 


Br(6xp) 


D+ - 




2.51 


3.08±0.36 


3.1 ±2.0 


o .1 m z.u 


q 07/1 _i_n f\Qa 


D° —, 


■ TT+K- 


4.03 


3.91±0.17 


5.8 ±0.7 


3.9 ± 1.0 


3.891±0.077 


D° - 


>Tr°K° 


1.35 


2.36±0.08 


2.1 ±0.6 


2.4 ±0.7 


2.38±0.09 


D° - 


-> K°r] 


0.80 


0.98±0.05 


0.9 ±0.2 


0.8 ±0.2 


0.96±0.06 


D° - 


*K°r,' 


1.51 


1.91±0.09 


0.3 ±0.2 


1.9 ±0.3 


1.90±0.11 


Dt- 


>K+K° 


4.79 


2.97±0.32 


5.1 ±0.9 


3.0 ±0.9 


2.98±0.08 


Dt 


-> iv + r) 


1.33 


1.82±0.32 


3.8 ±0.4 


1.9 ±0.5 


1.84±0.15 


D + s~ 


TT + rj' 


5.89 


3.82±0.36 


2.9 ±0.6 


4.6 ±0.6 


3.95±0.34 


D + s~ 















<0.06 



D — > PV decays, the intermediate states are pseudoscalar mesons with relatively large decay constants shown in Table 
HI The corresponding effective strong coupling constants are g p7rn = 4.2 obtained from p + — > 7r + 7r° or p° —> tt + it~, 
9K'Kit = 4.6 from K* (892)° — s- tt + K~, and g^KK — 4.5 from <f> — > K + K~ . For decays involving r] or r\\ we assume 
that g q — 4.2, g s = 4.6, and g ss = 4.5, where g q couple to the states with only u or d quarks, g s to two of the states 
with s quark, and g ss to all the three mesons with s quark, so that some effects of SU(3) breaking are considered. 

B. D — > PP 

For D —> PP decays, the decay rate is 

r(D^PP) = -^\A\ 2 , (34) 
87rm|j 

where p is the momentum of cither meson in the final state in the center-of-mass frame, p = 
\/(rri 2 D - (mp 1 + mp 2 ) 2 )(m z D - (m Pl -m^) 2 )/2m D . 

As is done in the naive factorization model, the Wilson coefficients a^'s are universal and process- independent, 
except with a relative strong phase for different topological diagrams. Because of the important nonperturbative 
effect of QCD in the charm system, a\ and ai should deviate a lot from the naive factorization approach. In order 
to give the most suitable results, we input the following values by hand, which also used the hint from the fit of 
diagrammatic approach: 

ai = 1.25 ±0.10, 

a 2 = (0.85±0.10) e '( 153±10 >°, 

a A = (0.90±0.10K( 160±10 >°, (35) 
a £ = (2.4±0.1)e l ( 55±10 >°, 

where a\ is the coefficient of the color-favored emission tree diagrams, ai for the color-suppressed emission diagrams, 
a a for the W-annihilation diagrams, and cle for the W-exchange diagrams. Large relative strong phases are considered 
in the a,, due to unneglected inelastic final-state interactions in the D decays. These values of a\ and are not far 
away from the large N c limit, except that we use quite large strong phases, which are required by the experimental 
data. As is discussed in the diagrammatic approach (i~9j . the W-annihilation contributions with hclicity-suppressed 
effect arc much smaller than those of the W-cxchangc diagrams. Therefore we use a much larger coefficient of cie than 
the W-annihilation coefficient a a- Besides, we ignore the disconnected hairpin diagrams, SE and SA, as discussed in 

The predicted branching ratios with annihilation-type contributions (all) and without annihilation type contribu- 
tions (emission) together with experimental measurements of charmed mesons decay into two pseudoscalar mesons 
are listed in Tables [TTTl IIV1 IVl for the Cabibbo-favored decays, the singly Cabibbo-suppressed decays, and the doubly 
Cabibbo-suppressed decays, respectively. There are many sources of theoretical uncertainties in the calculations. 
Since the decay constants of pseudoscalar and vector mesons are taken from experiments with very small errors, our 
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TABLE IV: Same as Tablelllll except for singly Cabibbo-suppressed decays of D — > PP(xlO 3 ). 



Modes 


Br(FSI) Br(diagrammatic) 


Br(emission) 


Br (all) 


Br(exp) 


D+ - 


+ TV+n 


1.7 


0.88±0.10 


1.0 ±0.5 


1.0 ±0.5 


1.18±0.07 


D+ - 


>K+K° 


8.6 


5.46±0.53 


11.3 ± 1.6 


8.4 ± 1.6 


6.12±0.22 


D+ 


-> Tv + r/ 


3.6 


1.48±0.26 


3.1 ± 1.0 


1.6 ±1.0 


3.54±0.21 


D+ - 


+ i 

-> 7T ?7 


7.9 


3.70±0.37 


3.7 ±0.7 


5.5 ±0.8 


4.68±0.29 


D° - 


S> 7T~ + 7T~ 


1.59 


2.24±0.10 


3.0 ±0.4 


2.2 ± 0.5 


1.45±0.05 


D° - 



-> 7T 7T 


1.16 


1.35±0.05 


0.7 ±0.2 


0.8 ±0.2 


0.81±0.05 


D° -> 


A+A" 


4.56 


1.92±0.08 


4.4 ±0.5 


3.0 ±0.8 


4.07±0.10 


D° - 


> K°K° 


0.93 








0.3 ±0.1 


0.64±0.08 


D° 


-¥ n°rj 


0.58 


0.75±0.02 


0.7 ±0.2 


1.1 ±0.3 


0.68±0.07 


D° - 




1.7 


0.74±0.02 


0.6 ±0.1 


0.6 ±0.2 


0.91±0.13 


D° 




1.0 


1.44±0.08 


1.3 ±0.4 


1.3 ±0.4 


1.67±0.18 


D° 


-> W 


2.2 


1.19±0.07 


0.04 ± 0.04 


1.1 ± 0.1 


1.05±0.26 


D + s~ 


>^°A+ 


1.6 


0.86±0.09 


0.9 ±0.2 


0.5 ±0.2 


0.62±0.23 


D + s~ 


>^+A° 


4.3 


2.73±0.26 


4.1 ±0.5 


2.8 ±0.6 


2.52±0.27 


Dt- 


-»■ A+7) 


2.7 


0.78±0.09 


0.8 ±0.5 


0.8 ±0.5 


1.76±0.36 




■> A+r/ 


5.2 


1.07±0.17 


0.7 ±0.3 


1.4 ±0.4 


1.8±0.5 



TABLE V: Same as Table ILUl except for doubly Cabibbo-suppressed decays of D — > PP(xW 



"Data from fH 
6 Data from [H 



Modes 


Br(diagrammatic) Br(emission) 


Br(all) 


Br(exp) 


D+ - 




1.98±0.22 


2.8 ±0.5 


1.7 ±0.5 




D+ - 


-> 7T K + 


1.59±0.15 


3.0 ±0.4 


2.2 ±0.4 


1.72±0.19 


D+ - 


K+T) 


0.98±0.04 


1.3 ±0.2 


1.2 ±0.2 


1.08±0.17 a 


D+ - 


-> K+T)' 


0.91±0.17 


0.4 ±0.1 


1.0 ±0.1 


1.76±0.22 6 


D° - 


■> 7T K 


0.67±0.02 


0.5 ±0.2 


0.6 ±0.2 




D° - 


> 7T-K+ 


1.12±0.05 


2.3 ±0.3 


1.6 ±0.4 


1.48±0.07 


D° - 


-»■ K°T) 


0.28±0.02 


0.23 ±0.05 


0.22 ±0.05 




D° - 




0.55±0.03 


0.08 ±0.06 


0.5 ±0.1 




D + s - 


> K+K° 


0.38±0.04 


0.7 ±0.4 


0.7 ±0.4 





numerical results are not very sensitive to the variations of meson decay constants. The branching ratios are truly 
sensitive to the coefficients of a\ and 02, especially to their relative strong phases. Since the systematic errors from 
theoretical models are usually difficult to estimate, we show uncertainties at the tables only from the parameters 
in Ea.([55|). For comparison we also show the recent results from the diagrammatic approach [l9[ and those consid- 
ering final-state interaction effects of nearby resonances (Hlj . It is clear that our results with large annihilation- type 
contributions agree with experiments much better than that of Ref . (Flj . For the Cabibbo-favored channels, which are 
the input data for \ 2 fit in the diagrammatic approach, we have comparable results with the diagrammatic approach 
[lfjl ]. For other channels, we have better agreement with experiments than the diagrammatic approach. The reason is 
mostly due to the SU(3) breaking effects, which had been fully neglected in the diagrammatic approach. 

The branching ratio of the pure annihilation process Df — > 7r + 7r° is vanished in our pole model. The resonant state 
that annihilates to 7r + 7r° is a scalar meson (0 ++ ), whose isospin could be 0, 1, or 2. However, isospin-0 would be ruled 
out because of charged final states, and isospin-2 is forbidden for the leading order AC = 1 weak decay. For the case 
of isospin-1, its G parity would be odd, which conflicts to a system of two pions whose G parity is even. Therefore, no 
resonant states can be produced and then annihilate to 7r + 7r°. In another word, no annihilation diagrams contribute 
to D + — > 7r + 7r and D + — s- 7r + 7r°. In fact, this kind of contribution is forbidden from the isospin symmetry of vt + 
and 7T° as identical particles. Simply, two pions can not form an s-wave isospin 1 state, because of the Bosc-Einstein 
statics. 

The pure annihilation process D° — > K°K°, with nonzero experimental data, also demonstrates the important 
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annihilation-type contributions. There are two kinds of contributions to this mode with dd and ss produced from 
weak vertex, respectively. In the flavor SU(3) limit, the rate vanishes due to the cancelation of CKM matrix elements, 
as predicted in the diagrammatic approach. Therefore, the effect of the SU(3) breaking is the dominant contribution 
here. In our pole model, we use /o(1370) and /o(1500) as two different poles with the dd and ss components, 
respectively, to describe the corresponding SU(3) breaking effect. We also refer to the argument of the long distance 
resonance effect in [l9|, [5f| , the t-channel final-state interaction in , the nonfactorizable chiral loop contributions 



57| . and the SU(3) breaking effect in the effective Wilson coefficients in |15| for this channel. 



Large branching ratios with 77' in the final states are both measured and predicted, which are larger than those 
with 77 in most cases, such as D° decays into K°r], K°r)' and into 77+77, 7r+7 ?'i although the phase spaces with 
77' are smaller than those with 77. For 77, the contributions from the components of dd and ss are destructive due to 
the minus sign in the mixing matrix of Eq. (|26[) and the positive mixing angle 1 ; while they are constructive for 77'. 
Besides, large W^-exchange contributions dominate most 77' modes, especially for D° — > ijrj' and D° — > K°r]', which 
demonstrates large annihilation-type contributions directly again. We also refer to this issue with i4Tg(1430) as an 
resonance in the spacelike form factors in the factorization ap pro ach in [52|, some effects of the inelastic final-state 
interactions in [58| , the final-state phases of the amplitudes in jl6l | , and the two-gluon anomaly effects in [59| . 



C. D -> PV decays 

The decay rate of D —> PV decays is 

T(D^PV) = ] JL_j:\A\>, (36) 

° pol. 

by summing over all the polarization states of the vector mesons. 

We assume that the coefficients cu are universal and process independent for D — > PV, but they are different from 
those of D — > PP as discussed in [3] and . Their absolute values are larger than those of D — > PP because the 
soft final-state interactions make more effects on D — > PV decays. In our calculations, they are used as 

af v = 1.32 ±0.10, 

ar = (0.75±0.10) e i ( 160±1 °)°, 

ar = (0.12±0.10)e*( 345±lor , (3?) 
a^ = (0.62±0.10)e i ( 238±1 °)°. 

Similarly to the PP modes, large relative strong phases due to inelastic final-state interactions are considered in 
the <2j. Again, the contributions from W-annihilation diagrams are smaller than the W-exchange ones. Besides, the 
relative strong phase between af v and a^ v is in accordance with the results from the diagrammatic approach [13 . [l9| . 

Our prediction of branching ratios of the Cabibbo-favorcd, the singly Cabibbo-suppresscd and the doubly Cabibbo- 
supprcsscd D —> PV decays are shown in Tables IVI1 1 VIII and IVIII1 respectively. The results in the third column 
(emission) in each of these tables are the predictions of rates with only the emission-type processes; while the results 
in the fourth column (all) also include the annihilation-type contributions. It is obvious that the annihilation-type 
contributions are of the same order as the emission- type diagrams, since the intermediate states here in the pole model 
are pseudoscalar mesons with relatively larger decay constants than those scalar mesons of the D — > PP case. Again, 
for theoretical uncertainty estimation, we use only those from the parameters shown in Eq.p7|), as illustration. For 
comparison, we also list the results of the diagrammatic approach [l9l | and the experimental date in these tables. It 
is easy to see that our results with the annihilation-type contributions agree with the experimental data. This means 
that the singlc-polc contribution dominates the annihilation- type contribution in most D — > PV decay channels. 
For example, although the D° — s- K°<p channel has no emission-type contribution, with vanishing branching ratio in 
the factorization approach, our pole model gives the right branching ratios agreeing with the experiment. This also 
confirms the calculation done in the perturbative QCD approach [201 ] . Besides, some of the SU(3) flavor symmetry 
breaking effects are considered in this work since the decay constants, transition form factors and effective strong 
coupling constants are involved. 

There is no resonant state contributing to the W-exchange diagram of D° — > ir°p° in the pole model, because a 7r° 
would violate the C parity, similarly to the case of D° — > rj{rj')uj, rj<p. Besides, the single-pole annihilation diagrams 



1 The theoretical and phenomenological estimates for mixing angle <f> is 42.2° and (39.3 ± 1.0)°, respectively |48H 
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TABLE VI: Branching ratios for Cabibbo- favored decays of D — > PV(%). The predicted rates with only emission-type 
contributions (emission) and with both annihilation- and emission-type contributions (all) are shown in the table, compared 
with the experimental data[2|J, the fitted results from the diagrammatic approach [li^] in which only the (A, Al) solution is 
quoted, and the results considering final-state interaction (FSI) effects of nearby resonances [5l|. 



Modes 


Br(FSI) Br(diagrammatic) Br(emission) 


JJI 1 dll J 


Br(sxp) 


u — 


> A p 


111Q 


i n Q4-9 o 
lU.ozrzZ.Z 


iz.z zn x.o 


O.O zt Z.Z 


lU.oztU. / 


D° - 


n. p 


fl RR 
u.oo 


1.54±1.15 


n 7 _i_ n c 

U.I I U.iJ 


1. / ± U. / 


-. qq+0.12 
i.oz_ 16 


D° - 


> n°K *° 


1 4Q 


o eo-i-n ^4 

Z . O Z; 1 1 J ■■ 


9 o _i_ 7 

Z.O ZL u. / 


9 o -I- i n 
z.y zt l.U 


z.ozztU.oO 


D° -t 


■ 7T + K*~ 


4.69 


91 +0 70 


3.8 ± 0.7 


o i i i n 
o.l =t l.U 


3- u °-0.53 


D° - 


-> r)K*° 


n ^i 




U.I zrz u. z 


n v i n 9 
U. f zt U.z 


U.yuztU.oU 


D° - 


* rj'K* 


U.UUiJ 


n m 9-1-0 flfH 


n 00^ -1- DDI 
u.uuo zn u.uui 


U.Ulu zt U.UUO 


<U.ll 


D° - 


-> K°uj 


2.16 


z . zuziz i.oo 


a + ^ 
U.U ZL u.u 


z.O zt U. / 


9 99-L-fl 1 9 

z.zzztU. Iz 


D° - 


H K°4> 




u.ouoiu.ioy 


o 

u 


U.o zt U.z 


U.oUoztU.UDU 


D+ - 




0.64 


1 83±0 49 


1.4 ± 1.3 


1 A 1 1 O 

1.4 zt l.o 


l.ODztU.lo 


D+ - 


+ K°p + 


11.77 


9.2±6.7 


15.1 ±3.8 


1 E 1 _l_ O Q 

±<j.± zrz o.o 


y .'izrzz .u 


D+H 


■ K+K*° 


3.86 




5.6 ± 1.9 


4.2 ± 1.7 


3.90±0.23 


D+H 


. k°K*+ 


3.37 




1.7 ±0.7 


1.0 ±0.6 


5.4±1.2 


Dt 


->■ W + 


9.49 




8.3 ± 1.3 


8.3 ± 1.3 


8.9±0.8[61J 


D + s 


-> Vp + 


2.61 




3.0 ±0.5 


3.0 ±0.5 


12.2±2.0 


Dt 


— > 


2.89 


4.38±0.35 


4.3 ±0.6 


4.3 ±0.6 


4.5±0.4 


D + s - 


-S-TT+p 


0.080 







0.4 ±0.4 


0.02±0.012 


D + s - 


^°p+ 


0.080 







0.4 ±0.4 




D + s 


— > 7r + aj 


0.0 










0.23±0.06 



can not contribute to the D — > pr], ttuj decays because of G parity violation. The isospin of resonant state for prj or 
7rw is one, so the G parity of the intermediate state is odd since it is a pseudoscalar meson. However, the G parity of 
p and T] are both even, and that of tt and u> are both odd, so the total G parity of the final states is even. Therefore, 
no resonant states are available for the decays of D mesons into pr\ and ttlu. It is even worse for the pure annihilation 
process Dj — > tt + uj, since its decay rate is predicted to be zero in the single-pole model, but it is not small in the 
experiment. This has already been discussed that this channel may be dominated by the final-state rescattering via 
quark exchange in (l9l. |55|. and by hidden strangeness final-state interactions in [60]. Besides, the pure annihilation 
mode Df — > tt + p° is predicted much larger in the pole model than the experiment data. The contributions from the 
two diagrams in this channel are constructive since the minus sign in the normalization of p° is compensated by the 
asymmetric space wave function of the two final states which are in the P-wave state. These two channels make such 
trouble that we fail to find a reasonable solution of Ap and Ay and predict the PV modes with the W-annihilation 
contributions in the diagrammatic approach fl9j. Hence, further discussions are still needed for these two channels. 

Dg — > pij is predicted much smaller than the mode of Dg — > prj } but the experimental branching ratios of the 
former is larger. This is a puzzle that the phase space of the former mode is much smaller than the latter, so its 
branching ratio should be smaller. In fact, the experimental measurement of Dg — > pif |(32j is already too old. It 
is already questioned by the PDG (26[, since this branching fraction (12.5 ± 2.2)% considerably exceeds the recent 
inclusive rj' fraction of (11.7 ± 1.8)%. 

IV. SUMMARY 

We have calculated the branching ratios for the two-body hadronic decays of charmed mesons into PP and PV 
using the generalized factorization approach for the emission-type diagrams and the pole-dominance model for the 
annihilation-type diagrams. Relative strong phases between different topological diagrams, which are important in 
the charmed decays, are considered in this work. Most of our predicted branching fractions are in accordance with 
the experimental data. Besides, compared to the naive and generalized factorization models ever before, the results 
in this work are much better since we have considered the annihilation-type diagrams and the relative strong phases 
between diagrams. 

We find that the annihilation-type contributions in the pole model are large for both PP and PV modes, which 
is also indicated by the difference between the life time of D + and Z?°. Comparing with the model- independent 
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TABLE VII: Same as Tablel VII except for singly Cabibbo-suppressed decays of D — > PV(x 10 3 ) 



Modes 


Br(FSI) Br(diagrammatic) 


Br(emission) 


Br(all) 


Br(exp) 


D° - 


■> 7T p + 


8.2 


8.34±1.69 


7.4 ±1.3 


10.2 ± 1.5 


9.8±0.4 


D°- 


■> 7T + p~ 


6.5 


3.92±0.46 


1.8 ±0.5 


3.5 ±0.6 


4.97±0.23 


D°- 


> „0 
-¥ 7T p 


1.7 


2.96±0.98 


1.4 ±0.6 


1.4 ±0.6 


3.73±0.22 


D° -»• 


K~K*+ 


4.5 


4.25±0.86 


5.5 ±0.8 


4.7 ±0.8 


4.38±0.21 


D° 


K+K*- 


2.8 


1.99±0.24 


2.0 ±0.3 


1.6 ±0.3 


1.56±0.12 


D° H 


■ K°K*° 


0.99 


0.29±0.22 





0.16 ±0.05 


<0.9 


£>° H 


■ K°K*° 


0.99 


0.29±0.22 





0.16 ±0.05 


<1.8 


D° 




0.08 


0.10±0.18 


0.08 ±0.02 


0.08 ±0.02 


<0.26 


D° 




1.1 


1.22±0.08 


1.0 ±0.3 


1.0 ±0.3 


0.76±0.05 


D° 




0.57 


0.31±0.10 


0.23 ±0.06 


0.23 ±0.06 


0.14±0.05 


D° 




0.24 


1.11±0.86 


0.05 ±0.01 


0.05 ±0.01 




D° 


-> r,'p° 


0.10 


0.14±0.02 


0.08 ±0.02 


0.08 ±0.02 




D° 


T)UJ 


1.9 


3.08±1.42 


1.2 ±0.3 


1.2 ±0.3 


2.21±0.23 


D° 


— » T)'uj 


0.001 


0.07±0.02 


0.0001 ± 0.0001 


0.0001 ± 0.0001 




D+ - 


-+1T+P 


1.7 




0.4 ±0.4 


0.8 ±0.7 


0.83±0.15 


D+ - 


->vT°p+ 


3.7 




5.3 ±1.7 


3.5 ± 1.6 




D+ —i 


• A+A*° 


2.5 




5.1 ± 1.1 


4.1 ± 1.0 


o 7f y+0.20 
°-' u -0.26 


D+ —i 


■ K°K*+ 


1.70 




14.0 ± 2.5 


12.4 ±2.4 


32±14 


D+ 


-J- r/p+ 


0.002 




0.4 ±0.4 


0.4 ±0.4 


< 7 


D+ 


' + 


1.3 




0.8 ±0.1 


0.8 ±0.1 


<5 


D+ 


-> 7T + 


5.0 


6.21±0.43 


5.1 ± 1.4 


5.1 ± 1.4 


5.44±0.26 


D+ 


->■ 7T + 0; 


0.35 




0.3 ±0.3 


0.3 ±0.3 


<0.34 


D + s - 


► 7T + A*° 


3.3 




2.3 ±0.8 


1.5 ±0.7 


2.25±0.39 


D + s - 


► 7T°A* + 


0.29 




0.4 ±0.2 


0.1 ±0.1 




Dt- 


+ K+p° 


2.4 




1.6 ±0.6 


1.0 ± 0.6 


2.7±0.5 


Dt- 


+ K°p+ 


19.5 




9.7 ±2.2 


7.5 ±2.1 








0.24 




1.0 ±0.4 


1.0 ±0.4 




Di- 


> r,'K*+ 


0.24 




0.4 ±0.2 


0.6 ±0.2 




D + s - 


-+ A+w 


0.72 




1.1 ±0.7 


1.8 ±0.7 


<2.4 




-> A+0 


0.15 




0.3 ±0.3 


0.3 ±0.3 


<0.6 



diagrammatic approach, wc reproduce their results with our specific model considering some SU(3) breaking effects. 
Furthermore, we get more predictions in many D — > PV decay channels, which are absent in the diagrammatic 
approach[l!|. Most of the results have a better agreement with experimental data than previous calculations. 
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TABLE VIII: Same as TablelVTl except for doubly Cabibbo-suppressed decays of D PV(xlO~ 4 ) 



Modes 


Br(diagrammatic) Br(emission) 


Br(all) 


Br(exp) 


D° H 


■ TY-K* + 


3.59±0.72 


3.7 ±0.6 


2.7 ±0.6 


o =4+1-80 


D° - 


> 7T K 


0.54±0.18 


0.6 ±0.2 


0.8 ±0.3 




D°- 




1.45±0.17 


1.1 ± 0.2 


0.9 ±0.3 




D° - 


^K°p° 


0.91±0.51 


0.2 ± 0.1 


0.5 ±0.2 




D° 




0.58±0.40 


0.2 ± 0.1 


0.7 ±0.2 




D° 




0.06±0.05 





0.20 ±0.06 




D° - 


■+ 7?A'*° 


0.33±0.08 


0.18 ±0.05 


0.17 ±0.05 




D°- 


4 r/A*° 


0.0040±0.0006 


0.001 ±0.001 


0.004 ± 0.001 




D+ - 


> TT+A"* 




3.0 ± 1.0 


2.2 ±0.9 


3.75+0.75 


D+ - 


» ^°A* + 




4.7 ±0.9 


4.0 ±0.9 




D+ - 


-S- K+p° 




1.4 ±0.4 


1.0 ±0.4 


2.1±0.5 


D+ - 


->K°p+ 




0.9 ± 0.4 


0.5 ±0.4 




D+ 


-> A+w 




1.4 ±0.5 


1.8 ± 0.5 




D+ 


-> K+cj) 







0.2 ±0.2 




D+ - 


-»■ ?7A* + 




1.5 ±0.2 


1.4 ±0.2 




D+ - 


4 r/A*+ 




0.013 ±0.006 


0.020 ± 0.07 




D + s~ 


► A+A*° 


0.20±0.05 


0.2 ± 0.2 


0.2 ±0.2 




D + s - 




1.17±0.86 


2.3 ±0.6 


2.3 ±0.6 





Appendix A: Individual formulas for various decay channels of D mesons 

The different contribution formulas for Cabibbo-favorcd decays of D — > PP are listed as 



A(D+ H- tt + K°) = i-j=V* s V ud [a 2 f K (m 2 D - ml)F^{m? K ) ± a x U(m 2 D - m 2 K )F DK (ml) 
A(D° -+ 7T+K-) = i%V: s V ud (arUiml - m 2 K )F^ K {m\) - a Egi f s f D ) , 

V2 \ m D ~ 771 Ki J 



A(D° -+ ir°K°) = i^fV c * s V ud [ a E9l f s f D "l D '""l ± a 2 f K (m 2 D - ml)F^{m 2 K ) 
2 \ m n — mi 



-.0 ,„o&Kh _ „• t/ -* t/ ( f _ f _ m D m K , f 2 

K 
,2 



A(D° -> = »%V£K d ( aa/jfCm?, - mg) ^ " g (m^ ) cos - asgifsfo ™ DmK \ (cos^-V^sin 



m D - m K , 

^4(D° -> A' ?/) = i-^V^Kd ( ^/^(mfj - m',)*^"' ( m D sin ^ ~ ^Egifsfo ^"'^ (sin + ^2 cos , 
2 V m D- m h 



A(D+ K+K°) = i^j=V* s V ud ( a 2 f K (m 2 Ds - m^)F D ^(m^) - a A9l f s f Ds J Ds ^ 2 



{ \_ TTL TJX \ 

-4(L>s -+7r + r/) = -iGjrV^* V ud ( -^=aif^{m 2 Ds - ml)F^ SVs (ml) sin + a A gif s fD s g Ds — ^cos^l, 
.A(£>+ -> tt+t/) = iG F F e *y Md f^ a i/^( m c s -mijF^^m^cos^-oxfli/s/pa ^ sma ° 2 sin A (Al) 

V V2 m D s - m a / 

where /g and g\ = 2.7 are, respectively, denoted as the decay constant of scalar mesons and effective strong coupling 
constant between the intermediate state and final states in the limit of SU(3) symmetry. Some phases from the 
propagators of the intermediate resonances are absorbed in the effective Wilson coefficients a a and a E - 
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The formulas for singly Cabibbo-supprcsscd decays of D — > PP are shown as 
A(D+^n+n°) = ^V^VMml-mDF^imDia.+a,), 
A(D + ^K+K°) = i ^(a 1 V; s V us f K ( m l-ml)F^ K ( m l)-a A V; d V ud9l fsf D 



-7= c 

A(D+ 7r+r?) = i^f cosTOKd^Cmi, - m^,f"(mj) - 2a A gifsl 
+V2a 2 [v; d V ud tf + V c * s V us f°](ml - ml)F°*(m 2 v )) , 



9 9 I ' 



m 2 r,m 



D"<-ao 



m D ~ m lo ' 



A(D+-Hr + rf) = i^(sm<f } V: d V u 4a 1 U(ml-ml,)F^(ml)-2a A g 1 f s f D 



AfnO . + -\ ;^ F T/*T/ /" / / 2 2\ v Dn/ 2\ if m D TO /o(1370) 

^ 7T 7r ) = «^=K cd Kd( aif v {m D - m n )F (roj - a E gifsjD 2 — g — 

m £ | m /o(1370) 



V2 

^^A ) = -&v: d vJa 2 Uml-ml)F^{ml) + a E9l f s f D ^"*f™>) \ 

1 \ m D m / (1370)/ 

A{D^K+K-) = l %V; s vJa 1 f K { m l- m ] < )F DK { m ] i )-a E g 1 f s f D ^f 50 ^ ) 

V2 V 771 D ~ TO /o(1500)/ 

</n0 , ^0^01 - Gf „ „ ^ f ( t/* t/ m J TO /o(1500) T/ * T ^ Wd to / (1370) 

.4(_L> ->• A A j = -l— -7=a,EgiJSjD\ VagVus— 2 ^~ cg£ *<iG0 2~ 



V2 V "» D -m^ riB001 m D -m fn 



D "7 (1500) "'£> "7o(1370). 

A(D° -> A) = i% fa 2 (m 2 D - m 2 )^™ 2 ^*^/,? + V^/'] - -^O^oa/^ - m 2 )^ (m 2 ) cos, 



,2 



-V2a E gifsfD 7 " £>ma ° cos0 



.Zip A') - (a 2 {m 2 D - m 2 )^(m 2 )[K* d K d /^ + K* s Ks/^] - ^=V c * d V ud a 2 U( m 2 D - m 2 )^* (m 2 ) sin0 

+V2a E g 1 f s fD " m 2 D ' maa sin 



^(I? ^) = J ^^ 2 (™ 2 D -m 2 )F o ^(m 2 )[y; d Kd/, d + ^K s / i ;]cos0 



TO D m /o(1500) . 2 , , 1 T ,* T/ TO l> m /o(1370) 2 



-aEgifsfD{V2V* s V us ^^ \ sin 2 <fi + —=V* d V ud — 2 



m D m 2 f n / 150( v, \/2 



cos 



-D "7o(1500) TO D TO / (1370) 



^(Z? -> W) = ^ (a 2 (m 2 D - m 2 )F^(m 2 n ,)[V c * d V ud f^ + V*V U8 f°,} coscf, 
+a 2 (m 2 3 - m 2 ,)F^( m 2 )[V^V ud f^ + V* s V us f*] sin^ 



-a Egifs f D (V2V; s V us ™ lmf f 500) sin20 - -^=V* d V ud f^f^ sin2(/) ) 



m / (1500) V 2 m -D TO / (1370) 

2 



^(D+^A+tt ) = -ffi (v^aAgifsfDs J Ds ™*l + K* d K rf a 2 /,(m 2 Js - m^F^K)) , 

z \ m Ds m K . j 

A(D+ -+ n+K°) = -& (v*V us a A gifsf Ds T^^f - ^K^i/^m 2 ^ - m^)^ 15 sK (ml)) , 

,2 



v 2 V to_d s - %■ V2 

-^(m 2 , - m|:)J^ x (m 2 v )[V; d V ud f^ + V c * s V us f°}\ , 

A(D+^K+r,') = (v*V us a Agi fsf Ds ^ mir f (J- sin^ + cos^) - y e * V us a x f K {m 2 D - m 2 v ,)F° srl °(m 2 K ) cos0 

V 2 V to £ , s - m^. V2 

-^(m 2 , - m^Jf^Cm^,)^^/^ + V C *K S /^]") , (A2) 
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The formulas for doubly Cabibbo-supprcssed decays of D — >• PP are listed as 



A(D+ -> A+tt 

^(L>+ -> 7T + A° 

A(D+ -> K+T) 

A(D+ -> K+T]' 
A(D° -> A+tt" 

^4(£)° -> AV 
„4(£>° -> A ?? 

A(D° -> A°r/ 
.A(£>+ -> A°A+ 



The formulas for Cabibbo-favored decays of £) — > PV are shown as 



^(A»° tt+A*- 
„4(A>° -> A"p+ 

^(A>° -> 7T°A* 

.A(£> -> A°p° 
^(A»° 77A* 
_4(L>° -> t?' A* 
_4(£>° -> A°w 
.4(L>° -> A°</> 

^(L>+ -> 7T + A*° 

_4(L>+ -> A'V 



„4(£>+ -> A+A 



.4(D+ -)■ A°A*+ 



-^KT d K s (lAgifsfD J ™*! + aif K (m 2 D - ml)F^{ml) 



'£) if,; 

,2 



-i%V£K. (aA9ifsfD ™ hmK ° - a 2 f K (m 2 D - ml)F°*(m 2 K ) 
y/Z \ m D - m K , 



G j? f TTXijyfTXj'C* i — Q 9 Dti 9 "\ 

i—V* d V us a A gifsfD— §— (V2sin<?!> - cos0) + a 1 f K (m D - m )F lq (m K ) coa</> , 

2 V m D- m KS ) 

™ D ™ K ° (sin0+ \/2cos</>) + a x f K {m% - m^)ff n ' (m 2 K ) sin A , 

l D m KX J 



G F 

t — V cd V us 



aAgifsfD- 



- l ^KdK s a E gifsfD— 5 aif K {m D - mjF (m K ) 

V 2 \ 13 - m A* 



c\ cd V US I L J d J U o 9 



^KdK,(a^i/ s /i3 1 



2 2 



-(•\/2sin0 - cos0) + a 2 /A-( m D - TO ^)-P(f' 7 ' ! ( m lr) cos^ . 



G ]? f 7TI J~)TTlp£* i— r. n Dl) 9 \ 

^-KrfKJ - a E gifsfD— §— (v2cos</> + sin0) + a 2 f K {m D - m n ,)F q (m K ) sin (j) \ , 

(A3) 



G 



= i^V* d V us f K (rn 2 Ds - m 2 K )F® sK {m 2 K ){ai +a 2 ). 



V2G F V c * s V u d(a E v g s f K fD 



' D 



+ a PV m K *f n A»"~(ml))(e*-p D ) 



DK* / 2 



L D _ A" 



V2G F V*V ud (a E v g s f K f D 



m 



D 



G F V* s V ud 
G E V* s V ud 
G F V* s V ud 
G E V* s V ud 
G F V*V ud 



'D "'A' 
2 

.PV „ m D 



+ a^m p f p F^{mi)\{e*-p D ) 



f V £ £ 

-a E g s f K f_ 



D ~ 2 2 



+ aZ v m K *f K *Ff*(m 2 K «)Ue* - PD ) 



PV t t 
— °E gsJKJD- 



D 



+ a^ v m p f K A^(m 2 K ))(e* ■ p D ), 



l D - K 



a E fKJD — J 



•d 



mr D -m' K 



(g s cos0- V2g ss sin0) + a% v m K * f K *F® Vq (m 2 K ,) cos0^ (e* -po), 



a E JKJD- 



' D 



- (g s sin 4> + V2g ss cos (f>) + a% v m K , f K , F^ Vq (m 2 K , ) sin <j> ) (e* ■ p D ) , 



l D K 



PV t f 
a E gsfKJD- 



V2G F V; s V ud a p E v g ss f K fD- 



m D -m K 



+ a^ muj f K A^(m 2 K ) (e*- PD ), 



■(£* -Pd), 



V2G F V c * s V ud m K , \a% v UA° K ' {ml) + a p 2 v ' f K *F^ : (m 2 K ,)j (s* ■ p D ), 
V2G F V c * 3 V ud m p (a PV f K A^(m 2 K ) + a PV f p F? K (m 2 p fj ( e * • p D ), 

+ a PV m K ,f K ,F 1 DsK (m 2 K ,))(e* -p Ds ), 



l Dc 



V2G F V*V ud a p A v g s f„f Ds — 



«A 9sfnfD s 



i> n 



2 9 ' "2 



mA*/A^ sK *(m^))(£*-Pi3 s ) 
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A(D+^r) P +) = -V2G F V^V ud arrn p f p F? sr, °(m 2 p )sm<t>(s* - PDs ), 
A(D+^ v 'p+) = V2G F V c * s V ud a PV m p f p F? s7 >°(m 2 p ) cos cf>(e* -p Ds ), 
A(D+ -> 7T+0) = V2G F V£V ud aF'm4UAg a *(ml)(e* -p Ds ), 

9 

771 

17)1 — 777 Z 



Ds 



TYl 

A(D+^ir°p+) = 2G F V c * s V ud a PV g q Uf Ds — , 

m D s - m i 



(e* - PDs), 



0, 



(A4) 



where the effective strong coupiing constants between the intermediate state and the finai states for D — > PV are 
g q = 4.2, g s = 4.6, and g ss = 4.5. The formulas for singly Cabibbo-suppressed decays of D — > PV are shown ~ 



as 



A(D° -»• TT+p- 

^(77° 7r°p° 
„4(L>° -> K-R*+ 
A(D° -> K+K*- 
A(D° -> A'°X*° 
„4(Z?° -> if A* 

.4(77° -> 7T CJ 

^(Z7° -> ttV 

.A(£> -> 77'p 
_4(.D -> 77W 
„4(L>° -> r/w 
A{D+ -> TT+p" 

„4(£>+ -> A >0 A* + 



— -G F V* d V ud I a 1 ^ 9qfnfD 



-V2a PV m p f p F^(m 2 ) )(e*-p D ) 



V: d V ud a PV m p [f p F^(m 2 p ) + f^ p (ml) ) (e* • pu), 



'7f 



\/2GfK* s K, ( g ss f D m 2 D {- 



J 'i 



V2G f V;* S K S ( a 1 ^ g ss f D m 2 D (- 



J 11 



f*> 



m D — ra~ m D — m~ 



) + a i m K* JkAq (m 2 K )\(e*-p D ) 



V2G F a pv f D m 2 D (v c * d V ud g s J* + V c 



s 



+ 



f s 

■I II 



rs V USySS 9 9 I 9 9 



V2G F a pv f D m 2 D [V c * d V u 



f* 



%d9s 9 9 



ss I 9 9 I 9 9 

m^, - mf m z D - mt 



(e* -Pd), 



(e* - Pd), 



G, 



V: d V ud a^ v m u (UF^(mi) - f v A»" (mi))(e* -p D ), 



= G F V^V us a% v m4 ,UF^(ml){e*-p D ) 

= G F V^V us a PV mtUF? 17 " (ml) cos 4>(s* -p D ), 

= G F a p 2 v m p ([V L * d V ud f* + V c * s V us f°}A£»(m 2 ) - l^V^V ud f p F^ (m 2 ) cos ^ (e* • p D ), 

= G F a PV m p ([V: d V ud tf, + V'V ua fy]A° p (m*, ) - -^V c * d V ud f p F^" (m 2 p ) sin A (e* ■ p D ) 

= G F a PV mu] ([V^f* + V c * s V U8 fs]A°»(m 2 ,) + -L cos <f>V^V ud f u F^ (m 2 )) (e* ■ p D ), 

= G F a PV m u {[V: d V ud ft r + V* s V us f^]A^ u (rn 2 ,) + -L sin cf>V c * d V ud f u F^ (ml )) (e* ■ p D ), 



G F V* d V ud m p [ 2a PV g q f^f D 



PV 



"D ""tt 



f p F?«(m 2 p ) a PV UA° p (ml) ) (e* ■ p D ) 



= G F V* d V ud m p [2a PV g q f v f D 



m n — m 



j a pv f p F^(m 2 p ) - a PV UA^(ml)) (s* -p D ), 

TT / 



V2G F ( V; d V ud a PV g s Uf D ™ D 2 + V c * s V us a pv m K «f K A° K ' (m 2 K ) (e* ■ p D ) 



m" D -m^ 



\f2G F [V; d Vuda P A v Qsf.fD- 



1 D 



V: s V us a pv m K *f K ,Fr K (m* K .) ) (e* ■ p D ) 



DK/2 



m n — m 
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A(D+ -> rjp + 
A(D+ -> 

.4(L>+ -> 7T+W 

A(D+ -> 7T°/v *+ 



.4(A>+ -> AV 



r/A 



*+ 



.4(13+ -> A+c 



A(D+ -> K~ 



A(D° -> 7T-1T+ 
.4(A>° -> 7r°A* 
A+/T 
^(D -> A /? 
A(D° -> A°w 
-4(D° -> tfV 

-4(£>° -> 7 ? A*° 
^(Z)° r/A*° 



= Gfiti p \V2a!- v A° p (ml){v; d V ud tf + V*V us f*} + V c * d V ud a p v f P F^ (m 2 p ) cos A (e* • p D ), 



= G F m p ( yfta™ A° p {m 2 v ,)[V: d V ud tf, + V c * s V as f*,} + V c * d V ad a p v f p F^" (m 2 p ) sin ^ ( e * • p D ), 

= V2G F V: s V us aP v m^^(ml){e* -p D ), 

= G F V* d V ud m u (a? v f u F?*(ml) + of v /^(m^ (e* • 

= vfe fa^o^./*/^ 2 mgs 2 + KdKdaf v m K >.UA° sK * (ml)) (e* ■ p Ds ), 
V mij. - mi, J 



G F [V: s V ua a p A v g s f K f Ds 



m 



m Ds - ™ 2 K 
2 

= G F ( V: s V us a p A v g s f K f Ds — ™ J 



V c * d V ud a PV m K ,UA° sK \ml)j (e* -p Ds ), 
V^V ud a PV m p f p F^ K (m 2 p )) (e* -p Dg ), 



V2G F (v; s V us a p A v g s f K f 



is ™ D _ S 2 + V; d V ud a pv m p f p Fr sK (m 2 p ) ) (e* ■ p Dg ), 

m D s K 



G F \ V L * s V U8 [a p A v f K f Ds 



m Ds ~ m K 



(g s cos 4> — V2g ss sin <fi) — \/2a p v uik* Ik* A , ] Ds ' ,s sin i 



+V2a PV m K .A° sK * (m 2 )[V* d V ud f* + V* s V us f*]j (e* ■ p Dg ), 

I 772 ^ 

g f V* s V us [a PV f K fD s ^2 — ° S i (9s sin </> + V2g ss cos<j)) + \pia PV m K ,j K *F® sr,s cos. 
m D s ~ m K 



-V2a PV m K *A° sK * (m 2 ,)[V* d V ud tf, + V*V us f°,} ) (e* ■ p Dg ), 



= ^2G F V c * s V us (a p A v g ss f K f 



IDs— 2 



G F ( V: s V us a PV gJ K f Ds f Ds 2 + V* d V ud a PV m u f u F? sK {m 2 J ) (e* ■ p Dg ). 



L D S 



+ a^m^UF^im 2 ) + a pv m ^f K A^(m 2 K )) (s* -p Ds ), 



(A5) 



The formulas for doubly Cabibbo-supprcsscd decays of D —> PV are shown as 



V2G F V c * d V us (a PV g s f K f D ™» 2 + a PV m K ,f K > F^{m 2 K .)) (e* ■ p D ) 



m D ~ m K 



= G F V^V us (a PV m K ,f K ,F^(m 2 K ,) - a PV g s f K f D 



1 D 



m D- m K 



(e* -Pd), 



V2G F V c * d V us f K (a pv m p A° p (m 2 K ) - a pv gJ D ™ D 2 ) (e* ■ p D ) 



G F V c * d V us f K (a PV m p A° p (m 2 K ) - a PV g s f 



in 



D 



D — 2 T 

m D — m K 

2 
D 



(e* - Pd), 



G F V; d V us f K (a PV m u A^{m 2 K ) - a pv g s f D ™» a ) (e* ■ p D ), 

\ m D m K J 



~2 

= V2G F V^V us a pv g s f K fD ™ D 2 (e* -Pd), 

m D — m K 



= GFV* d V us a F v f K fD 



G F V* d V us a F v fxfo 



1 d - " l K 



(g s cos0- V2g ss s'mcj)) + a PV m K * j K *F^ q im 2 K ») cos^ (e* ■ p D ), 
(g s sin<f) + V2g ss cos./)) + a PV m K ,j K , F® Vq (m 2 K ,) sin 0^ (e* -p D ), 
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A(D+ -> 7T+A'* ) 
^(L>+ -> 7r°A*+) 
.A(£> + -> 

-> a+^) 

^(L>+ -> AT+0) 
.4(L>+ -> »/#*+) 
^(L>+ -> 7/AT* + ) 



V2G F V* d V us (a p A V 9 s f K fD 2 m ° a + afm^./K.^lm^)) (e* ■ pr>), 
G F K* d K s (a p /g a f K fD 2 m2 ° 2 - a PV m K ,f K ,F^{m 2 K ,)\ (s* ■ p D ), 

\ m D ~ m K J 

GFV* d V us f K ( a PV g s fD 2 — : - 2 - a i Vm p A o P (' m K) ) ( e * ■ Pd), 



f Z5 

\/2Gf V* d V us ftc (a PV g s f D ™ u 2 + af^^K)) (e* ■ Pu), 
GFV&VnfJaFg.fD 2 m ° 2 + a PV m aJ A^(m 2 K ))(s* -p D ), 



2 



(2 
a P A !k!d 2 _° 2 (9s cose/) - \/2g ss sin0) + a pv m K *f K *F[ 
m D m K 

GFV* d V us I a PV f K f D mD (g s sine/) + V2g ss cos0) + a pv mic* f k* F® nq {m 2 K ,) 



„4(Z?+ K+K*°) = V2G F V c * d V us m K * (a 2 f K .Ff> sK \m 2 K .) + a pv f K A° sK * (m 2 K ^j (s*-p D ), 
V2G F V; d V us m K , (a PV f K *F? sK {m 2 K .) + a PV f K A^ K * {m 2 K )\ (e* ■ p D ). 



A(D+ K°K*+) 



(m^,) c os0|(e* -pr>), 
$ I (e* • Pd), 

(A6) 
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